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ABSTRACT 

Time  spreading  measurements  provide  an  indirect  measure  of  the  acoustic 

bandwidth  tha^can  be  supported  by  the  water  channel,  which  is  critical  to  the  design  of 

sonar  sv  ?ems  Thne  spreading  measurements  were  collected  in  a  water  channel  100  m 
sonar  systems,  lime  sprc  s  Data  wgre  collected  at  frequencies  of  20-22  kHz,  27- 

^kHz  and  35-37  kHz  using  linear  FM  pulses  2  s  in  duration.  The  experiments  were 
part  of’  a" 'collaborative  TOP  trial  known  as  Trial  Scotiar ^ 

Source-receiver  and  target  position  were  recorded  using  a  portable  tracking  rang 
a  u  tko  t tq  Tn  the  mner  time  spreading  measurements  are  compared  with 
modelled  Stim  teTob.  usTn^eteSna?  Mode,  (GSM).  The  GSM  estimares 

7 time SSm  due  to  multipath  propagation  compare  favourably  with  the 
experimental  data§  However,  time  spreading  of  individual  paths  -  bey°nf*®  > 
predictions  -  is  also  evident.  Additionally,  the  data  were  used  to  compute  the  standard 
deviation  of  the  received  echo  intensity  at  each  frequency.  The  standard  deviation  was 
rnmnuted  two  different  ways.  First  it  was  computed  using  the  peak  echo  level  from  each 
KrSfteLpv  Then,  it  was  computed  from  the  total  energy  received 
from  each  of  te  pings.  At’all  three  frequencies,  the  standard  deviation  was  lower  by  1  to 
2  dB  when  computed  from  the  total  received  energy. 

RESUME 

Les  mesures  d’etalement  dans  le  temps  fournissent  une  mesure  indirecte  de  la 
larveur  de  bande  acoustique  que  pent  accepter  le  canal  de  son  dans  l'eau  ce  qm  est  une 
dnnnee  critique  pour  la  conception  des  systemes  sonar.  Les  mesures  d  etalement  dans 
temps  ont  et6  recueillies  dans  un  canal  d’une  profondeur  de  100  m au riarge  de  la  cote  de 
Nnnvelle-Ecosse  Les  donnees  ont  ete  recueillies  aux  frequences  20  22  kH  , 
km  et  35  37  kHz,  en  utilisant  des  impulsions  FM  lineaires  d’une  duree  de  2  s.  Les 

experiences  faisifnt  partie  d’un  essai  de  collaboration  du  programme  de  cooperation 
technique  TTCP  appele  Trial  Scotian  (HF)  (essai  HF  de  la  Nouvelle-Ecosse),  orgaiiis 
par  le  groupe  de’spdcialistes  de  l’acoustique  environnementale  du  P™Srarranetechmque 

GTP-11,  contre-mesures  et  armes  sous-marines.  Le  Canada,  sonar  actif 

Rnvaume-Uni  ont  participe  a  cet  essai.  Le  reseau  canadien  SEAHORSE,  sonar  actii 
haute  frequence  ancre,  servait  de  source  et  de  recepteur,  tandis  qu  un  repeteur  d  echo  a 
derive  l^breT servait  de  cible.  Les  positions  de  la  source  et  du  recepteur,  d  une  part,  et 
cedes  de  la  cible  d’ autre  part,  ont  ete  enregistrees  a  l’aide  d  une  installation  de  poursuite 
Soloitle  par  les  e“u.  La  communicalion  compare  les  mesures  d-etalement  aux 
Smalions  ob  enues  au  moyen  du  logiciel  Generic  Sonar  Model  -  GSM  (modele  de 
sotT Srique)  fpartir  des  donnees  gdomdtriques  de  l’expdrience.  Les  esumations 
GSM  de  Semen?  dans  le  temps  dfia  la  propagation  sur  divers  trajets.se .compare 
favorablement  aux  donnees  experimentales  Toutefois,  pour  des  trajets  indmduels,  on 
constate  aussi  un  etalement  dans  le  temps  depassant  celut  qut  etatt  preyu  par  le  GSM  De 
n?ns  tes  donndes  ont  servi  it  calculer  L&art  type  de  l’intensite  de  1’echo  re{U  a  chaque 

freauence.  L'ecart  type  a  ete  calcule  de  deux  manteres  dtfferentes,  d  abord  a  partir  du 

frenuence  donnee,  puis  a  partir  de 


i  partir 
1’ autre  methode. 
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EXECUTIVE  SUMMARY 


Background 

It  is  well  known  that  increasing  the  frequency  bandwidth  of  a  sonar  pulse  can 
increase  the  spatial  resolution,  and  therefore  improve  target  localization.  One  might 
mistakenly  argue  then,  that  increasing  the  sonar's  bandwidth  will  always  enhance  its 
performance.  However,  the  response  of  the  water  channel  to  the  sonar  signal  places  an 
upper  limit  on  the  improvement  that  one  may  achieve.  In  effect,  the  water  channel  itself 
is  bandwidth  limited. 

Time  spreading  measurements  provide  an  indirect  measure  of  the  acoustic 
bandwidth  that  can  be  supported  by  the  water  channel,  which  is  critical  to  the  design  of  a 
sonar  system.  The  sonar  designer  uses  this  information  to  optimize  the  transmit  pulse 
employed  in  the  sonar.  Decisions  such  as  the  pulse  length,  the  pulse  type,  and  the 
frequency  bandwidth  are  all  based  in  part  on  the  time  spreading  characteristics  of  the 
water  channel. 


Results 

Time  spreading  measurements  were  collected  in  a  water  channel  100  m  deep,  off 
the  coast  of  Nova  Scotia.  Data  were  collected  at  frequencies  of  20-22  kHz,  27-29  kHz, 
and  35-37  kHz  using  linear  FM  pulses  2  s  in  duration.  The  experiments  were  part  of  a 
collaborative  TTCP  trial  known  as  Trial  Scotian  (HF)  organized  by  the  Environmental 
Acoustic  specialists  group  of  GTP-11,  Underwater  Weapons  and  Countermeasures. 
Canada,  the  US,  and  the  UK  participated  in  the  trial.  The  trial’s  location  reflects  the 
increased  interest  in  littoral  waters  expressed  by  the  Canadian  Navy. 

In  the  paper,  time  spreading  measurements  are  compared  with  modelled  estimates 
obtained  using  the  Generic  Sonar  Model  (GSM).  The  GSM  estimates  of  time  spreading 
due  to  multipath  propagation  compare  favourably  with  the  experimental  data.  However, 
time  spreading  of  individual  paths  -  beyond  the  model  predictions  -  is  also  evident.  The 
results  of  this  study  will  be  used  to  improve  performance  predictions  of  present  day 
sonars  as  well  as  to  aid  in  the  design  of  future  sonars. 
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INTRODUCTION 


This  paper  presents  measurements  of  time  spreading  over  the  frequency  band  20 
kHz  -  40  kHz,  collected  in  a  shallow  water  channel  off  the  coast  of  Nova  Scotia.  Time 
spreading  measurements  provide  an  indirect  measure  of  the  (water)  channel  bandwidth 
which  is  critical  to  the  design  of  sonar  systems.  The  measurements  are  compared  with 
predictions  obtained  using  the  Generic  Sonar  Model 1  (GSM). 

The  experimental  measurements  were  part  of  a  collaborative  TTCP  trial  known  as 
Trial  Scotian  (HF)  organized  by  the  Environmental  Acoustic  specialists  group  of  GTP- 
11,  Underwater  Weapons  and  Countermeasures.  Canada,  the  US,  and  the  UK 
participated  in  the  trial.  Additional  measurements  made  during  the  trial  from  DREA's 
high-frequency  SEAHORSE  array2  included  submarine  target  strength,  surface  ship 
target  strength,  high  frequency  ambient  noise,  shallow  water  reverberation,  and  frequency 
spreading.  The  results  from  these  experiments  as  well  as  (water)  channel  characterization 
measurements  made  with  a  US  General  Test  Vehicle  (GTV)  will  be  reported  separately. 

Following  the  introduction  we  describe  the  experiment  and  the  bottom  topography 
of  the  experimental  site.  Then  the  time  spreading  data  are  presented  and  compared  with 
estimates  obtained  using  GSM.  Finally,  we  present  some  conclusions  which  resulted 
from  the  research. 

I.  EXPERIMENTAL  EQUIPMENT  AND  GEOMETRY 

Figure  1  shows  a  schematic  of  the  geometry  for  the  time  spreading  experiment. 
The  UK  echo  repeater  was  deployed  to  76  m  depth;  SEAHORSE  was  deployed  to  42  m. 
The  seabed  was  nominally  flat,  with  an  average  water  depth  over  the  extent  of  the 
experiment  of  104  m.  The  benthic  topography  is  discussed  in  greater  detail  in  the 
following  section.  The  experiment  was  conducted  inside  the  perimeter  of  a  portable 
tracking  range  developed  and  operated  by  personnel  from  Naval  Undersea  Warfare 
Center  (NUWC),  Division  Keyport.  SEAHORSE  was  decoupled  from  surface  motion 
using  its  two-stage  decoupling  system2  and  anchored  to  the  seabed  within  the  tracking 
range  using  a  second  two-stage  decoupling  system.  The  experiment  was  conducted 
during  a  2  hour  period  throughout  which  the  echo  repeater  was  free  drifting.  This 
allowed  for  time  spreading  measurements  at  several  SEAHORSE-to-target  ranges. 
Figure  2  shows  a  plot  of  the  SEAHORSE  and  the  echo  repeater  GPS  tracks  measured 
during  the  experiment.  The  start  locations  of  the  tracks  are  as  marked.  Figure  3  contains 
a  plot  of  relative  range  and  bearing  from  SEAHORSE  to  the  echo  repeater  over  the 
duration  of  the  experiment.  The  data  in  Figure  3  equate  to  an  average  drift  speed  for  the 
echo  repeater  of  0.1  m/s. 

A.  The  SEAHORSE  System 

SEAHORSE  is  a  HF  active  sonar  for  collecting  environmental  acoustic  data  in  the 
open  ocean.  The  sonar  operates  over  the  frequency  band  20  kHz  to  40  kHz  and  has  an 
approximately  symmetrical  beamwidth  of  0.1  steradians.  In  conjunction  with  the 
acoustic  sensors,  SEAHORSE  is  instrumented  with  a  range  of  non-acoustic  sensors  to 
assist  in  the  evaluation  of  the  data.  The  non-acoustic  sensors  include  depth,  tilt,  roll,  and 
heading  sensors,  to  monitor  array  position  and  direction,  as  well  as  accelerometers  to 
monitor  platform  motion.  A  series  of  motors  contained  in  the  sonar  head  -  and  remotely 
controlled  from  a  research  ship  -  point  the  sonar  head  to  the  required  azimuth  and  tilt. 
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Figure  1:  Schematic  of  experimental  setup.  Note  that  the  echo-repeater  was  free  drifting 
during  the  experiment  thereby  varying  the  range  from  1300  m  to  2000  m. 
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Figure  2:  Tracks  of  SEAHORSE  and  echo  repeater  during  the  experiment. 
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Figure  3:  Relative  range  and  bearing  from  SEAHORSE  to  the  echo 
experiment. 


repeater  during  the 


n!nUla4:  i  Contour  Plot  °f the  water  depth  (in  meters)  at  the  experimental  site.  The  dot 

^^^JfTEAHORSF^  rf]adingsfrom  whif  the  contours  were  extrapolatea 
l  lie  location  oj  SEAHORSE  during  the  experiment  is  denoted  by  an  X. 


Communication  between  the  research  ship  and  SEAHORSE  is  handled  via  two  RF  links. 
The  first  is  a  one-way  UHF  data  link  operating  at  2.37  GHz  which  is  used  to  tiansmit  the 
acoustic  and  non-acoustic  data  collected  by  SEAHORSE  back  to  the  ship.  The  second  is 
a  two-way  VHF  command  link  operating  at  462  MHz  which  is  used  to  maintain  system 
control  (e.g.  select  pulse  type,  set  amplifier  gains,  control  pointing  direction,  e  c.). 

The  transmitter  consists  of  a  48-element  array  of  transducers,  shaded  such  that 
one  obtains  an  approximately  conical  beam  of  0.1  steradtans,  “V^sured  to 
nnintQ  fin  cross-section  this  equates  to  a  two-sided  beamwidth  of  20  .)  Beam  steering 
s  accomplished  through  purely  mechanical  means.  The  sonar 
freauencv  band  20  kHz  to  40  kHz  with  a  source  strength  of  202  ±  3  dB  re  4tFa  <s>  lm. 
The  same  48-element  array  used  for  the  acoustic  transmitter  is  used  as  the  primary 
a^usb^recdve^As^receiver,  the  48-element  array  of  transducers  is  shaded  to  obtain 
an  approximately  conical  beam  of  0.1  steradians,  measured  to  the  -3  dB  points^  The 
receive  sensitivity  is  -180  ±  5  dB  re  1  pPa/Hz  across  the  frequency  band  20  to  40  kHz. 


B.  The  DEMONSTRATOR  Echo  Repeater 


DEMONSTRATOR  is  a  free-drifting  echo  repeater  designed  and  operated  by  the 
UK's  Defence  Research  Agency  (DRA)  System  control  £ 

command  link  operating  at  (nominally)  400  MHz.  DEMONSTRATOR  can  be 
programmed  to  operate  at  center  frequencies  of  22,  28,  or  37  kHz  with  a  bandwidth  at  any 
of  these  center  frequencies  of  4  kHz.  The  target  strength  is  also  programmable  and  was 
set  to  +10  dB  during  the  entire  experiment.  Both  the  transmittei  and  the  iecei\e 
hydrophone  are  omnidirectional.  As  well  as  returning  a  duplicate  of  the  received  signal 
DEMONSTRATOR  returns  a  19  m/s  Doppler  shifted  version  of  the  signal  following  a 
delav  of  200  ms.  The  system  contains  an  onboard  DAT  recorder  with  a  tape  storage 
capacity  of  two  hours.  The  system  has  a  battery  life  of  4  to  6  hours  continuous,  but  can 
be  shut  down  and  reactivated  as  desired  to  extend  system  life. 

II.  ENVIRONMENTAL  DESCRIPTION 

A.  Benthic  Topography 

The  site  of  the  experiment  was  the  eastern  flank  of  the  La  Have  Bank,  about  100 
nautical  miles  south  of  Halifax,  Nova  Scotia.  Figure  4  contains  a  contour  plot ^of  the 
bottom  for  the  experimental  site.  The  dots  overlaid  on  the  plot  mark 
depth  measurements  taken  using  the  bathymetnc  system  on  CF  Q  ,  Hpnth 

ship  employed  for  the  sea  trial.  The  contour  lines  were  constructed  fiom  these  depth 
measurements.  The  location  of  SEAHORSE  is  denoted  in  the  figure  by  an  X.  The 
bathymetric  measurements  were  also  used  to  generate  the  surface  plot  of  the  benthic 
interface  presented  in  Figure  5.  These  plots  show  that  the  bottom  had  a  gentle  slope  to 
the  southeast,  with  the  depth  increasing  from  100  to  1 10  m  over  about  10  km  Although 
there  is  little  small-scale  variability  in  the  bathymetry,  the  bottom  at  the  site  is  rough 
Side  scan  sonar  and  visual  records  collected  at  the  site  during  a  subsequent  trial  showed 
the  bottom  at  the  site  to  be  rough  to  very  rough.  The  bottom  is  predominately  cobble, 
with  the  smoothest  areas  consisting  of  course  gravel.  In  the  roughest  areas,  cobbles  weie 
the  order  of  0.3  to  3.0  m  diameter. 

B.  The  Water  Column 

Figure  6  shows  the  sound-speed  profile  for  the  site  taken  shortly  after  the 
experiment.  The  minimum  speed  occurred  at  42  m  depth,  coinciding  (appioximately) 
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Figure  5:  Surface  plot  of  the  seabed  at  the  experimental  site.  The  figure  was  constructed 
from  the  echosounder  readings  shown  in  Figure  4. 


speed  of  sound  (m/s) 

Figure  6:  Sound-speed  profile  of  the  water  column  taken  shortly  after  the  experiment. 


Figure  7 •  The  10  dominant  eigenrays  between  the  SEAHORSE  transmitter  and  the  echo 
repeater  for  the  21  kHz  geometry.  For  clarity,  rays  departing  the  source  toward  the 
surface  are  solid  lines  and  rays  departing  the  source  toward  the  bottom  are  dashed.  The 
bold  solid  line  is  the  most  intense  ray  and  the  bold  dashed  line  is  the  first  ai  rival. 


range  (m) 


Figure  8:  Target  returns  from  the  echo  repeater  for  two  "typical"  linear  FM  pings 
compared  to  the  idealized  output  of  the  matched  filter  (solid  line)  illustrating  time 
spreading  in  the  received  echoes. 
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with  the  SEAHORSE  depth.  Figure  7  contains  the  eigenrays  from  SEAHORSF  to  the 
echo  repeater  for  the  10  most  intense  rays,  for  a  source  frequency  o”  2 1  kHz  In  mis  a^d 
subsequent  figures,  a  nominal  sound  speed  of  1470  m/s  was  used  to  convert  iime  un  tsm 
Ifrnfted  As  f  ne  rf"  See  fr°'u  lhe  ray  Pallls.  the  sound-speed  profile  of  Figure  6  resulted  in 
geometry  “  b0U"Ce  Pa‘hS  fr°m  SEAH0RSE  ®  echo?epea,er  for^ex^rimemS 

III.  THE  TIME  SPREADING  EXPERIMENT 

™  LSSSStL  \  L™d  S 

paper,  each  pulse  type  is  referred  to  by  its  center  frequency:  21,  28,  or  36  kHz. 

ss!FFrVF““ 

between  the  start  of  the  transmit  pulse  and  reception  of  the  echo  reneLr  t  n™  c  Y 
dropouts  m  the  data  telemetry  link  -  due  to  higPh  seas  -  resuhed  TatlsTs^  LFM 

so  tha/th?**  Stquenre:  Analysis  was  performed  on  sixteen  pulses  from  each  sequence 
amh?  nui?ber  of  senes  examined  was  consistent  at  all  frequencies.  Durino 
y?.is  PuIses  were  heterodyned  down  in  frequency  to  the  band  500  9500  h7  i  a 
pass  filtered,  and  sub-sampled  by  a  factor  of  8  to  redVce  •  250b  Hz’ low' 

samples1  (5 0%  ^  verlapf be^wee^s^mples)1^ ^d ^FourieMransfomecl 

SM2 

0  m/s.  It  may  appear  somewhat  odd  to  convert  time  spreading  data  to  ranee  units 

A.  Time  Spreading  at  21  kHz 

the  21  nifl8  contains  a  plot  of  the  time  spreading  returns  from  two  LFM  pulses  for 
the  21  kHz  run.  Each  return  is  normalized  to  0  dB  at  its  respective  maximum  The  two 

P?k  rCtUrn  fr0m  the  e"h° 

Swnin  te^f 1400  "  ThC  °Pehn  Hne 

that  the  data  are  time  spread  significantly  below  the  -4  dB  level.  § 
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range  (m) 


Figure  9:  Time  spreading  data  for 
target  returns,  measured  relative  to 


the  maximum,  mean  (see  footnote  5),  and  minimum 
the  dominant  eigenray  path,  for  each  of  the  16  pulses 


at  21  kHz. 


range  (m) 


Figure  10: 
relative  to 
at  21  kHz. 


Time  spreading  data  for  the  mean  target  return  (see  footnote  5),  measured 
he  total  energy  returned  by  the  target  via  all  paths,  for  each  of  the  16  pulses 
For  companion  the  mean  level  measured  relative  to  the  dominant  eigenray 


path  is  replotted  from  Figure  9. 
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To  indicate  the  variability  in  the  received  echoes.  Figure  9  shows  a  plot  of  the 
time  spreading  data  for  target  returns  from  three  of  the  21  kHz  LFM  pulses.  The  data 
plotted  are  the  echo  levels  for  the  maximum,  the  minimum,  and  the  mean5  target  returns 
as  measured  relative  to  the  dominant  eigenray  path  (i.e.  the  peak  level  of  the  curve)  for 
the  21  kHz  pulses.  Using  the  peak  level  from  each  of  the  21  kHz  LFM  pulses  one 
obtains  a  standard  deviation6  of  2.3  dB.  Alternatively,  one  may  describe  the  strength  of 
an  echo  in  terms  of  the  total  energy1  returned  by  the  target  via  all  paths.  Using  this 
definition  one  obtains  a  standard  deviation  of  1.7  dB  for  the  16  pulses  at  21  kHz.  Figure 
It)  contains  a  plot  of  the  data  for  the  mean  target  return  as  measured  relative  to  the  total 
energy.  For  comparison  the  mean  level  measured  relative  to  the  dominant  eigenray  path 

Lr9Pi  ?«  * ?/i!frUr1 9'  TablC  1  HstS  the  Standard  deviation  using  both  approaches  for 

the  2 1 ,  28,  and  36  kHz  data  sets. 

TABLE  I.  Standard  deviation  (a)  of  the  target  echo  for  linear  FM  (LFM)  pulses  in 
individual  frequency  bands.  The  standard  deviations  are  computed  for  echo 
strength  expressed  as  peak  amplitude  and  total  received  energy. 


Frequency  (kHz) 

2l 

28 

36 


a  (dB)  using  peak  level 
23 

6.5 

4.5 


G  (dB)  using  total  energy 

TT 

4.7 

2.2 


Figure  1 1  contains  a  plot  of  the  average  of  the  time  spreading  data,  normalized  to 

^e'0d?  at I thC  2/  kHz  run'  Prior  t0  waging,  the  matched-filtered  data 
were  time  shifted  to  align  the  peak  return  from  the  echo  repeater  for  each  ping.  The  onen 

ScC  ,es  repres^t  the  afrage  of  all  16  pings.  The  solid  line  in  the  figure  shows  the  spatial 
esolutmn  of  the  matched  filter.  This  curve  was  generated  by  using  an  exact  duplicate  of 
the  kernel  as  the  input  signal.  To  remain  consistent  with  the  data,  the  solid  line  was 
obtained  by  averaging  the  matched-filter  output  of  16  realizations  of  the  input  signal 
each  input  signal  having  a  slightly  different  time  delay.  In  this  way,  any  smearing  of  the 
data  resulting  from  misalignment  of  the  FFT  window8  and  the  received  echoes  will  be 
included  in  the  ideal  matched-filter  curve.  The  averaged  data  substantiates  the 
observation  made  earlier  in  reference  to  the  individual  pings  show  in  Figure  8.  That  is: 
little  or  no  time  spreading  occurs  in  the  region  of  the  peak,  down  to  about  -4  dB  but 

KJTT.Sf  time  SPread‘ng  increases  substantially  Furthermore,  the  averaged  data 
indicates  that  the  spieading  is  asymmetric  about  the  maximum. 

B.  Time  Spreading  at  28  kHz 

of  the  S?taiI!S  a  pl?.°f  the  tim®  Heading  data  for  target  returns  from  three 

of  the  28  kHz  LFM  pulses.  The  range  from  SEAHORSE  to  the  echo  repeater  was 

approximately  1860  m.  The  data  plotted  in  the  figure  are  the  matched-filter  outputs  for 
the  maximum,  the  minimum,  and  the  mean  echoes,  as  measured  relative  to  the  dominant 
eigenray  path,  for  the  28  kHz  pulses.  (Recall  footnotes  4  and  5.)  It  is  obvious  from 

l2h°  FigUre  9  that  there  is  §reater  variability  in  the  received  echoes  at 
fs  ii  ’r  S  21  kf^'  ln,  ™ore  quantitative  terms,  one  obtains  a  standard  deviation  of 
,  A  dB  fiom  the  peak  levels  from  the  sixteen  28  kHz  pulses  compared  to  2  3  dB  at  21 
kHz.  Using  the  total  energy  returned  by  the  target  via  all  paths,  one  obtains  a  standard 
eviation  of  4.7  dB  for  the  16  pulses  at  28  kHz.  Figure  13  contains  a  plot  of  the  28  kHz 
data  for  the  mean  target  return  as  measured  relative  to  the  total  energy  (The  curve 
representing  the  mean  level  measured  relative  to  the  dominant  eigenray  path  is  not 
reproduced  from  Figure  12  because  at  28  kHz  it  happens  to  be  the  sam!  ping  ) 
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kHz.  run. 


at  28  kHz- 
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F*gwre  13:  Time  spreading  data  for  the  mean  target  return  (see  footnote  5)  measured 
at  28  kHz.  €  mer8y  retUmed  by  the  target  Vm  allPaths’ for  each  of  the  16  pulses 


- ideal  filter 

(16  simulations  averaged) 
°  data  level 

(16  pings  averaged) 


28  kHz 


1  -10 


XT  1  \Ll  ■  1  ' 

1850  1855  I860  1865  1870  1875 

v.  ..  ,  range  (m) 

tigure  14:  Averaged  time  spreading  data  for  16  pings  ( open  circles)  normnlb *4  tn  n  An 
28mZnmm'  C0mpared  ,0  the  ldeallzed  OU,pm  0f,he  mal cbed  filter  (solmZe)fcfhe 


Figure  14  contains  a  plot  of  the  averaged  time  spreading  data,  normalized  to  0  dB 
^  r-  tup  o q  kHz  run  The  solid  line  in  the  figure  is  the  output  of  the 

„h  ?pnfater  for  each  ping  The  open  circles  represent  the  average  of  all  16  pings. 

asymmetric  about  the  maximum  and  the  most  significant  time  spreading  occurs  before, 
not  after,  the  dominant  eigenray  arrival. 

C.  Time  Spreading  at  36  kHz 

Fienre  15  illustrates  the  variability  of  the  time  spreading  measured  at  36  kHz. 

The  figure  contains  the  matched-filter  outputs  (with  peaks  aligned)  for  the  maximum 

•  ‘mnm  mrl  mean  target  returns  as  measured  relative  to  the  dominant  eigenray  path  fo 
minimum  and  mean  targ  jhe  f  SEAHORSE  to  the  echo  repeater  was 

three  of  the  36  each  of  the  36  kHz  LFM  pulses,  one 

approximately  •  §  TJsin^  the  total  energy  returned  by  the  target  via 

eigenray  path  is  replotted  from  Figure  15. 

Figure  17  contains  a  plot  of  the  averaged  time  spreading  data,  normalized  to  zero 
&  n  o/t  \r\jrj  rnn  Xhp  solid  line  in  the  figure  is  the  output  of  the 

-8  dB  level. 

In  the  following  section  we  compare  the  averaged  time  spreading  data  with  results 
modelled  using  the  Generic  Sonar  Motfe!  (GSM)  to  gain  some  tnstgh,  tnto  the  cause  of 
the  measured  time  spreading. 

IV.  COMPARISON  OF  TIME  SPREADING  DATA  TO  GSM 

Time  spreading  for  the  experimental  geometry  was  modeM  Mg  GSM.  Inputs 

reneater  When  running  the  model  we  assumed  a  pulse  length  of  0.67  ms  to  oDtai 
tenlporai  resofmion  equivalent  to  that  of  the  LFM  pulses  employed  in  the  experiment. 

T,  v-inttnm  hackscatterin^  strength  was  estimated  using  a  Lambert's  rule 
depending  wT  a  27  dB^  independent  of  f„_  Reverbemtm" 

lack  of  Quantitative  bottom  information  currently  available  for  the  site  * 
containsqthe^ bottom  loss  coefficients  employed.  These  numbers  represent  reasonable 
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valueS  and  provide  3n  acceptable  fit  to  the  data.  Note  that  the  bottom  reflection 
coefficient  at  grazing  angles  greeter  th3n  8°  does  not  affect  the  results  for  this  geometry. 

TABLE  II.  Bottom  reflection  coefficient  table  employed  for  GSM  evaluation.  Note  that 

the  bottom  reflection  coefficient  at  grazing  angles  greater  than  8°  does  not 
ajject  the  results  for  this  geometry. 


_ _ Grazing  angle  (degrees) _ _ Bottom  reflection  coefficient  (dB) 


uu  -0.00 

2-°  -2.00 

4-°  -4.00 

60  -6.00 

8-°  -8.00 

100  -8.00 

20-°  -8.00 

900  -8.00 

MTTT  The  GSM  predictions  were  obtained  by  calculating  the  dominant  eigenrays  (using 
MULTIP)  from  SEAHORSE  to  the  target,  (in  this  case  the  echo  repeater)  to  obtain  a 
senes  of  arrivals  at  the  target  for  a  specific  SEAHORSE-to-target  range.  The  source  level 
and  arrival  time  of  each  of  these  eigenrays  were  then  used  as  sources1 1  at  the  echo 
lepeater  and  propagated  back  to  SEAHORSE.  The  signal  excess  subroutine  CMP  AX  1  in 

aAHnDcr°df  „m  order  to  output  the  echo  level  plus  reverberation  level  received  at 
SEAHORSE.  Finally,  the  eigenray  arrival  times  were  converted  to  range  using  a  sound 
speed  of  1470  m/s  so  as  to  remain  consistent  with  the  data.  The  ranges  of  the  eigenray 

e^enrays1  (recall  Rgure  7)  ^  th0USht  °f  HS  approximate  path  len§ths  for  the  various 

As  indicated  in  Figure  3  the  echo  repeater  opened  from  SEAHORSE  at  a  rate  of 
approximately  0.1  m/s  during  the  experiment.  This  drift  rate  resulted  in  a  change  of 
range  of  30  to  80  m  (1%  to  3%  of  the  total  range)  over  the  time  taken  to  transmit  the  pin° 
sequence  at  each  frequency.  The  sensitivity  of  the  modelled  output  to  these  range 
vanatmns  is  indreated  m  Figure  18  in  which  the  GSM  results  at  28  kHz  for  ranges  of 
1860  m,  1875  m,  1890  m,  and  1905  m  are  shown.  The  four  curves  were  time  shifted  to 
align  the  maximum  of  each  curve.  The  ambient  noise  level  lies  at  approximately  -25  dB 
Clearly,  the  amplitude  and  number  of  the  peaks  varies  with  even  small  changes  in  range 
If  the  modelled  returns  are  averaged,  as  was  done  with  the  experimental  data,  the  peaks 

Swre  m  m  Smgie  ™dellec!  returns  wil1  be  smeared.  To  include  the  range  averaging 
effects  in  the  modelled  results,  16  GSM  estimates  were  taken  at  each  frequency  to 

range^fnrV^W  16  expenmental  Pin§s-  For  simplicity,  we  calculated  the  start  and  stop 
range  for  each  frequency  sequence  from  the  measured  time  delays  and  assumed  a 

dunnThrahfhP^  1°  hIrTh  ^  intermediate  ran§es-  The  measured  increase  in  range 
2  d  28  kHz  exPenEents  was  30  m  and  thus  the  model  estimates  were 
mputed  at  2  m  range  increments.  The  measured  increase  in  range  during  the  36  kHz 
experiment  was  75  m  and  thus  the  model  estimates  were  computed  every  5  m. 

Recall  from  Table  I  that  the  standard  deviation  computed  from  the  peak 

anSil  th65  m  th?  5  P1,nf  varied  f!om  2'3  t0  6-5  dB-  These  values  are  considerably 
g  eate1  than  the  standaid  deviation  of  the  16  model  estimates,  which  was  approximately  1 

exne  if  SUenC1?-  T,hat  1S  t0  S3y’  ping-to-ping  intensity  fluctuations  measured 
experimentally  are  largely  unaccounted  for  in  GSM.  Therefore,  when  computing  the 

average  time  spread  from  the  data,  each  ping  was  normalized  to  zero  dB  at  its  maximum, 
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range  (m) 


Figure  15:  Time  spreading  data  for  the  maximum,  mean  (see  footnote  5),  and  minimum 
target  returns,  measured  relative  to  the  dominant  eigenray  path,  for  each  of  the  16  pulses 

at  36  kHz ■ 


range  (m) 


Figure  16:  Time  spreading  data  for  the  mean  target  return  (see  footnote  5),  measured 
relative  to  the  total  energy  returned  by  the  target  via  all  paths,  for  each  of  the  16  pulses 
at  36  kHz.  For  comparison  the  mean  level  measured  relative  to  the  dominant  eigenray 

path  is  replotted  from  Figure  15. 
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- ideal  filter 

(16  simulations  averaged) 

o  data  level 

(16  pings  averaged) 


range  (m) 

Figure  17:  Averaged  time  spreading  data  for  16  pings  ( open  circles),  normalized  to  0  dB 
at  its  maximum,  compared  to  the  idealized  output  of  the  matched  filter  ( solid  line),  for  the 
36  kHz  run. 
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Figure  18:  GSM  results  at  28  kHz.  for  ranges  of  1860  m,  1875  m,  1890  m,  1905  m.  The 
four  curves  were  time  shifted  to  align  the  maximum  of  each  curve. 


prior  to  averaging.  This  meant  that  each  ping  had  an  equal  weighting  and  ensured  that  no 
single  ping  dominated  the  average,  thereby  providing  a  more  effective  compaiison  with 
GSM.12  In  the  subsequent  discussion,  we  shall  refer  to  the  data  normalized  to  zero  dB  at 
its  maximum  prior  to  averaging  as  the  weighted  data,  and  the  data  normalized  to  zero  dB 
at  its  maximum  after  averaging  as  the  unweighted  data.  Note  that  the  data  shown 
previously  in  Figures  11,  14,  and  17  were  unweighted. 

A.  GSM/Data  Time  Spreading  Comparisons  at  21  kHz 

Figure  19  contains  a  comparison  of  the  average  of  the  modelled  echoes  for  the  21 
kHz  sequence  with  both  the  weighted  and  unweighted  average  for  the  21  kHz  data,  for  a 
nominal  SEAHORSE  to  echo  repeater  separation  of  1400  m.  Squares  represent  the 
weighted-average  data,  open  circles  represent  the  unweighted-average  data,  and  the  solid 
line  represents  the  model  estimate.  The  GSM  eigenray  output  (not  shown)  indicates  that 
the  direct  path  is  neither  the  dominant  arrival  nor  the  earliest  arrival.  It  is  rather 
labourious  (and  not  particularly  instructive)  to  associate  each  modelled  peak  with  its 
respective  eigenray  due  to  the  combined  effect  of  the  two-way  propagation  employed  in 
the  model  and  the  averaging  of  the  model  output.  However,  it  is  worth  mentioning  by 
way  of  example  that  the  dominant  eigenray  for  the  21  kHz  data  is  the  ray  that  departs 
SEAHORSE  upward  (i.e.  toward  the  surface),  then  refracts  downward  and  finally  refracts 
upward  before  arriving  at  the  echo  repeater.  This  eigenray  is  the  emboldened  solid  line  in 
Figure  7.  For  comparison,  the  first  arrival  is  the  single  bottom  reflection  path  shown  as 
the  emboldened  dashed  line  in  Figure  7. 

Clearly,  there  is  less  than  perfect  agreement  between  the  model  and  the 
experimental  data.  The  model  does  however,  explain  the  dominant  feature  in  the  data 
—  the  asymmetry  of  the  data  about  its  peak  amplitude.  This  feature  is  consistent  with  the 
distribution  of  peaks  predicted  by  GSM.  That  is  to  say,  GSM  predicts  far  fewer  arrivals 
prior  to  the  dominant  peak  than  after.  Although  secondary  peaks  are  not  visible  in  the 
data  the  relative  amplitudes  of  the  peaks  in  the  model  are  consistent  with  the  general 
decay  in  amplitude  of  the  data  after  the  arrival  of  the  dominant  peak.  Furthermore, 
recalling  Figure  7,  the  dominant  peak  in  Figure  19  has  no  surface  or  bottom  interaction 
whereas  the  peak  at  1395  m  as  well  as  peaks  located  at  1405  m  and  beyond  have  bottom 
and/or  surface  interactions.  This  provides  phenomenological  evidence  at  least  that  one 
might  alter  the  bottom  reflection  coefficients  in  Table  II  to  reduce  these  secondary  peaks 
without  reducing  the  dominant  one,  thereby  improving  the  fit  between  GSM  and  the  data. 
However,  without  quantitative  information  on  the  bottom  parameters,  the  exercise  would 
not  be  particularly  instructive.  More  importantly,  any  improvement  in  the  fit  is  unlikely 
to  account  for  the  discrepancy  between  the  data  and  the  model  in  the  immediate  vicinity 
of  the  dominant  peak.  That  is  to  say,  there  is  significant  time  spreading  of  the  dominant 
data  peak  beyond  that  attributable  to  multipath  propagation.  It  seems  unlikely  that  the 
spreading  in  the  data  results  from  ocean  surface  motion  since  the  model  predicts  no 
surface  interaction  for  the  dominant  path.  Nor  is  platform  motion  likely  to  be  responsible 
since  both  SEAHORSE  and  the  echo  repeater  were  nearly  motionless. 

B.  GSM/Data  Time  Spreading  Comparisons  at  28  kHz 

Figure  20  contains  a  comparison  of  the  average  of  the  modelled  echoes  at  28  kHz 
with  both  the  weighted  and  unweighted  average  for  the  28  kHz  data.  The  nominal 
separation  from  SEAHORSE  to  the  echo  repeater  during  this  ping  sequence  was  1860  m. 
Squares  represent  the  weighted  average  of  the  data,  open  circles  represent  the  unweighted 
average  of  the  data,  and  a  solid  line  represents  the  model  estimate.  As  in  the  21  kHz  case, 
the  structure  of  the  data  is  reasonably  well  modelled  by  GSM.  The  small  undulations  in 
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range  (m) 

Figure  19:  Averaged  time  spreading  data  for  the  16  pings  at  21  kHz,  compared  to  the 
averaged  output  of  the  Generic  Sonar  Model  for  16  ranges  at  2  m  increments  ( solid  line). 
Both  a  weighted  data  average  (squares),  computed  by  normalizing  to  0  dB  prior  to 
averaging,  and  an  unweighted  data  average  (circles),  computed  by  normalizing  after 
averaging  are  shown.  (Refer  to  text.) 

the  data  from  1845  m  to  1853  m  may  correspond  to  the  three  small  peaks  in  the  model  at 
1846,  1849,  and  1852  m.  Both  data  and  model  are  significantly  more  symmetric  about 
the  maximum  than  at  21  kHz.  The  28  kHz  data  are  time  spread  relative,  to  the  ideal 
matched  filter  (recall  Figure  14).  However,  the  model  results  suggests  that  multipath 
propagation  should  have  resulted  in  even  greater  time  spreading  than  that  exhibited  by  the 
data.  However,  this  discrepancy  between  GSM  and  the  data  may  result  from  inaccuracies 
in  the  bottom  loss  and  backscattering  strength  tables  as  discussed  previously. 

C.  GSM/Data  Time  Spreading  Comparisons  at  36  kHz 

Figure  21  contains  a  comparison  of  the  average  of  the  modelled  echoes  at  36  kHz 
with  both  the  weighted  and  unweighted  data,  for  a  nominal  SEAHORSE  to  echo  repeater 
separation  of  1950  m.  Squares  represent  the  weighted  average  of  the  data,  open  circles 
represent  the  unweighted  average  of  the  data,  and  a  solid  line  represents  the  model 
estimate.  Once  again  the  structure  of  the  data  is  reasonably  well  modelled  by  GSM.  The 
peaks  in  the  model  that  occur  from  1935  m  to  1940  m  roughly  correspond  to  peaks  in  the 
evenly  weighted  data.  Note  that  none  of  these  peaks  are  clearly  visible  in  the  unweighted 
data.  As  with  the  28  kHz  case,  the  model  predicts  more  extensive  time  spreading  about 
the  main  peak  than  is  exhibited  by  the  data. 

V.  DISCUSSION 

The  energy  calculation  of  the  standard  deviation  indicates  that  there  is  significant 
ping-to-ping  variability  in  the  energy.  This  likely  occurs  as  ping-to-ping  fluctuations 
redistribute  the  energy  into  paths  which  don’t  intercept  the  target,  or  as  more  energy  is 
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Figure  20:  Averaged  time  spreading  data  for  the  16  pings  at  Z6  ktLz>  compareaxo  me 
overused  output  of  the  Generic  Sonar  Model  for  16  ranges  at  2  m  increments  (solid  line  . 
Both  a  weighted  data  average  (squares),  computed  by  normalizing  to  0  dB  prior  to 
averaging,  and  an  unweighted  data  average  (circles),  computed  by  normalizing  after 
averaging  are  shown.  (Refer  to  text.) 

scattered  incoherently  through  surface  and/or  bottom  interactions  from  one  ping  than 
from  another  No  attempt  was  made  to  model  this  because  GSM  lacks  a  realistic  fading 
model  for  individual  eigenrays.  That  is  to  say,  the  GSM  results  do  not  have  any 
statistical  fluctuations  built  into  the  eigenray  calculation  -  partly  because  GSM  is  not  a 
range-dependent  model. 

At  all  three  frequencies,  the  standard  deviation  was  lower  when  computed  using 
the  total  energy  than  when  computed  using  the  peak  value.  (Recall  Table  I.)  This  is 
because  the  peak  value  calculation  results  from  a  single  eigenray  and  small  changes  in 
ran^e  affect  the  amount  of  energy  in  the  dominant  eigenray.  In  contrast,  the  energy 
calculation  is  less  sensitive  to  these  fluctuations  since  a  significant  portion  of  the  energy 
stripped  from  one  eigenray  will  reach  the  receiver  via  a  slightly  different  path(s).  Note 
that  the  difference  in  the  value  of  the  standard  deviation  obtained  from  the  energy 
calculation  and  the  peak  level  calculation  is  in  good  agreement  with  the  standard 
deviation  of  1  dB  obtained  from  the  model,  which  was  cited  m  the  previous  section. 

Although  GSM  provides  a  qualitative  explanation  of  the  time  spreading,  there  are 
definite  discrepancies  between  the  model  and  the  data.  For  example,  the  21  kHz  data  in 
Figure  19  exhibits  significantly  more  time  spreading  immediately  after  the  arrival  ol  the 
dominant  peak  than  estimated  by  the  model.  Also,  in  Figures  20  and  21  the  several  well 
defined  peaks  around  the  maximum  obtained  from  GSM  are  not  all  observable  in  the 
data  It  is  reasonable  to  assume  that  the  agreement  would  improve  if  one  had  more 
accurate  inputs  into  GSM  for  the  bottom  loss,  the  backscattenng  strength,  and  the 
SEAHORSE-to-target  separation  for  each  ping.  However,  the  range-independent  nature 
of  GSM  will  clearly  limit  one's  ability  to  match  the  fine  structure  of  the  data. 
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Figure  21:  Averaged  time  spreading  data  for  the  16  pings  at  36  kHz,  compared  to  the 

ZT8  0UtPufflhe  Genenc  Sonar  Model  for  16  ranges  at  5  m  increments  (solid  line) 
Both  a  weighted  data  average  (squares),  computed  by  normalizing  to  0  dB  prior  to 

avei  aging,  and  an  unweighted  data  average  (circles),  computed  by  normalizing  after 
averaging  are  shown.  (Refer  to  text.)  normalizing  after 

VI.  SUMMARY  AND  CONCLUSIONS 

r  KT  Time  spreading  measurements  were  collected  in  water  100  m  deep,  off  the  coast 
of  Nova  Scotia.  Data  were  collected  at  center  frequencies  of  21  kHz,  28  kHz  and  36  kHz 

SEAHORSE  FM  PU  S6S  V61  ln  duration  with  a  2-kHz  bandwidth!  Canada's 
SEAHORSE  array,  an  anchored,  high  frequency  active  sonar  was  emploved  for  the 

source-receiver,  and  a  UK  free-drifting  echo  repeater  was  employ*  f£?the  target 

b^hCeeirrtre^rSet  P0Sltl°n  were  rec;orded  usin§  a  portable  target  range  operated 
.v/i?  US-  At  each  frequency,  time-spreading  measurements  obtained  from  averaging 

ModeHGSMT  AHH rgS  Trefrmpar,ed  estimates  obtained  the  Generic  Sona? 

*e  pi®s  was lomprf’  Sta"dard  de™ti0n  °f  the  received  echo  in""sity  from 

Firct  -t The  standard  deviati°n  of  the  received  echo  intensity  was  computed  two  wavs 
list  it  was  computed  from  the  peak  echo  level  from  each  of  the  16  pulses  at  a  given 
frequency  Then,  it  was  computed  from  the  total  energy  received  from  each  of  the  16 
pngs'  ^  a  thl?e  frequencies,  the  standard  deviation  was  lower  by  1  to  2  dB  when 
' d  fromth®  total  received  energy.  This  is  likely  because  the  maximum  value 
results  from  a  single  eigenray  and  the  energy  in  the  dominant  eigenray  closely  depends 
n  small  changes  in  range,  etc.  In  contrast,  the  energy  calculation  is  less  sensitive  to 

slmhtlyTfferenTpShr  3  Significant  P°rtion  of  the  ener§y  win  sti11  reach  the  receiver  via 
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At  each  freauencv  the  average  time  spreading  was  computed  and  normalized  to 
0  dB  with  to  the  maxtaum.  These  data  were  compared  to  the  output  of  an  .deal 

matched  filterwith  the  following  results.  At  21  and  28  kHz,  little  or  no  time  spreading 
"  in  the  region  of  the  PV  down  .o  about  the  4  dB  po.n^  Mow ,4 1  dB 
however,  the  time  spreading  increases  substantially.  At  36  kHz,  the  time  spreading  on  y 
becomes  significant  below  the  -8  dB  points. 

GSM  provides  an  acceptable  qualitative  explanation  of  the  time  spreading, 
narticularlv  in  describing  the  overall  width  and  shape  of  the  mam  returns.  At  21  kHz,  the 
lime  spSi“g  data  wert  extremely  asymmetric  about  the  maximum ,m .the  < lata  — 
at  98  and  36  kHz  the  data  were  somewhat  more  symmetric.  The  GSM  results  in 
Lfthe  degrefof  symmetry  resulted  from  the  arrival  time  and  amplitude  of  the  various 
eigenravs  §There  were  some  definite  discrepancies  between  the  model  and  the  data  At 
21  ml  the  data  exhibited  significantly  more  time  spreading  in  the  vicinity  of  he 
nraGmnm  than  estimated  by  the  model.  At  28,  and  36  kHz  the  GSM  output  exhibited 

several  clear  peaks  which  were  not  well  defined  in  the  data.  It [he 
agreement  to  improve  if  one  had  more  accurate  inputs  into  GSM  for  the  bottom  loss,  me 
EsSng  strength,  and  the  SEAHORSE-to-target  separation  for  each  pmg. 
However,  the  range-independent  nature  of  GSM  will  clearly  limit  one  s  ability  to  match 
the  fine  structure  of  the  data  . 
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